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INTRODUCTION 


\  .  ..  I 

^  Our  work  during  this  period  continued  to  have  as  its  main 

objective  the  study  of  (a)  the  thermodynamics  and  structure  of  plasmas 
and  dense  molecular  fluids  and  (b)  the  distribution  of  the  microfield 
in  a  strongly  coupled  plasma.  Our  research  included  both  rigorous 
analysis  and  approximation  schemes. 

We  continued  applying  our  medium  approximation  to  a  variety 
of  systems.  This  method  continues  to  be  the  simplest  scheme  available 
for  dealing  with  dense  diatomic  fluids.  It  also  seems  to  work  reason¬ 
ably  well  for  simple  mixtures.  In  fact,  we  found  a  new  scheme  for  re¬ 
ducing  a  mixture  to  a  one  component  fluid,  which  we  hope  to  explore 
further. 

In  the  area  of  plasmas  our  APEX  aoproximation  seems  to  he 
the  best  one  around  for  static  microfield  distribution.  We  have. 


during  the  last  year,  obtained  some  new  results  about  the  time  depen¬ 
dent  gas,  but  no  simple  approximation  scheme  yet.  _ _ _ 

Our  work  has  been  in  close  collaboration  with  groups  at 
Los  Alamos  National  Laboratory  and  at  Lawrence  Livermore  National 
Laboratory.  We  give  here  a  brief  description  of  the  work  and  attach 
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1.  On  Potential  and  Flold  Fluctuations  In  Clnaslcnl  CharKod 
Systems 

(J.  L.  Lobowltz  and  Ph.  A.  Martin) 

Using  clcccroscacic  Idcncicics  cite  potential  and  microfield  in  a 
plasma,  important  for  determining  line  sitapes,  are  expressed  as 
limits  of  local  quarTtities.  These  are  shown  to  be  well  defined 
for  typical  configurations  of  macroscopic,  i.e.  infinite  systems 
(under  some  mild  clustering  assumptions).  Their  covariance  con¬ 
tain  a  slowly  decaying  part  (  |XI”^  .  the  potential)  whose 

coefficient  is  universal  whenever  the  S t  i  I  I  ingc r-Lovc 1 1  second 
moment  condition  holds.  U'c  show  further  that  the  contributions 
from  distant  regions  (which  arc  equal  to  suitable  averages  over 
local  regions)  have  a  Gaussian  distribution. 

(Appeared  in  Jour.  Stat.  Phys. ,  Vol.  34,  Nos.  1/2,  1984)  Reprints 
under  separate  cover 

2.  Low-Frequency  Component  Electric  Microflold  Distributions  in 
Plasmas 

(Carlos  A.  Iglesias,  Joel  L.  Lebowitz,  David  MacGowan  and  Hugh 
E.  DeWitt) 

We  evaluate  the  low-frequency  component  electric 
microfield  distribution  at  a  charged  point,  i.e.  the  field 
distribution  in  a  gas  of  ions  interacting  through  a  Yukawa 
potential.  The  method  employed  is  an  adaptation  of  the  APEX 
approximation  previously  developed  for  the  high-frequency 
component  and  involves  a  nonintcracting  quasi-particle  represen¬ 
tation  of  the  screened  plasma  designed  to  yield  the  correct  second 
moment  of  tho  mlcrofleld  distribution.  The  APEX  results  are  compared 
to  new  Monte  Carlo  simulations  at  various  plasma  coupling  strongth.s 
and  we  find  good  agreement.  In  connection  with  this  analysis  we 
also  solvod  the  hypernattsd  chain  approximation  for  tho  Yukawa 
potential  and  compared  that  to  the  Monte  Carlo  elmulatlona. 

(To  appear  in  Phye.  Rev.  A,  teat.  pub.  date  March  1985)  Reorinte 
under  separate  cover  when  availab'*^ 
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3.  Thornodyaamics  of  Komonucloar  Diatomic  Fluids  from  tho  Angular 
Median  Potential 

(David  MacGowan,  Eduardo  M.  Waisman,  Joel  L.  Lobowitz  and 
Jerome  K.  Percus) 

The  use  of  the  angular  median 

pocenclal  a  Cemparacure-lndependanc  spharlcal  rataranca  ayacam  lot  ap- 
proxLmaClng  Mlacular  fluids  la  CaaCad  for  lea  pradlcclona  of  charmodyna- 
■Ics.  Calculaclons  hava  baan  carrlad  out  for  a  wlda  ranga  of  hooionuclaar 
dlaComlcs  ulch  conclnuous  acoa-acom  pocanclala  ballavad  Co  ba  raprasancaclva 
of  cha  full  ranga  of  sloHilaclon  data  avallabla  for  such  syscaas.  Tha 
results  for  cha  prassura  ara  surprisingly  good  boch  in  cha  daconaclon  regime 
and  around  cha  crlpla  point.  In  cha  laccer  case,  however.  Che  Incernal 
energies  for  highly  alongacad  molecules  wlch  accracctva  pocenclal  wells  are 
considerably  too  positive.  Comparison  with  other  perturbation  theories 
indicates  that  the  median  reference  system  gives  better  pressures  but 
poorer  energies  chan  RAM.  and  chat  In  aany  cases,  especially  for  purely 
repulsive  potentials.  It  gives  results  of  comparable  accuracy  to  chose 
obtained  with  non-spher leal  reference  systems. 

(Appeared  in  J,  Chem.  Phys.  80(0),  15  March  1984)  Reprints  under 
separate  cover 

4.  Spherical  Reference  Systems  for  Kon-Spherical  Hard  Interactions 
(G.  0.  irilliains,  J.  L.  Lebowitz  and  J.  K.  Percus) 

We  investigate  the  applicability  of  the  median  and  Darker- 
Henderson  prescriptions  for  obtaining  spherical  reference  systems 
for  three  models:  hard  linear  triatomlca,  hard  heteronuolear  duihb- 
bells,  and  tup-component  mlxturen  of  hard  duiid>b«lls.  We  propose 
na  empirloal  method  for  determiainf  the  median  potential  for 
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systems  lacking  a  high  degree  of  symmetry.  For  mixtures  of  hard 
molecules,  we  find  that  both  the  median  and  Barker-Hendorson  pros¬ 
criptions  give  rise  t.o  approximately  additive  hard  sphere  roforonco 
potentials. 

(Appeared  in  J.  Phys.  Chero. ,  December  1984)  Reprints  under  separate 
cover  when  available 

5.  Monte  Carlo  Simulation  of  Hard  Spheroids 

(J.  W.  Perram,  M.  S.  fferthelm,  J.  L.  Lebowitx  and  G.  0.  Wllllnms) 
We  present  Monte  Carlo  simulations  of  the  equation  of  state 
and  radial  distribution  function  for  a  model  fluid  composed  of  linrd 
spheroids. 

(Appeared  in  Chem.  Phys.  I.ctts.  ,  Vol.  105,  No.  3,  16  March  1084) 
Reprints  under  separate  cover 

6.  Electric  Microfield  Distributions  in  Multicomponent  Plasmas 
(Carlos  A.  Iglesias  and  Joel  L.  Lebowltz) 

We  evaluate  the  electric  microfield  distribution  in  a  multi- 
component  plasma  (MCP).  The  method  employed  is  an  adaption  of  the 
very  simple  adjustable-parameter  exponential  approximation  previ¬ 
ously  developed  for  one-component  plasmas  (OCP).  We  also  discuss 
a  still  simpler  approximation  in  which  the  MCP  is  replaced  by  an 
effective  OCP.  The  results  are  generally  close  to  each  other  and 
the  former  Is  in  very  good  agreement  with  computer  simulations. 

(Appeared  In  Phys.  Rev.  A,  Vol.  30,  No.  4,  October  1984)  Reprints 
under  separate  cover 
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7.  Ono  Molecular  Fluid  Approximation  for  Diatomic  Fluid  Mlxturoa 
(Eduardo  M.  Waisman,  Joel  L.  Lobowitz  and  David  MacGowan) 

Ue  Investigate  a  one  component  molecular  fluid  approximation  for 

conformally  similar  molecules.  Ue  test  this  scheme  on  (two)  mixtures 

of  rigid  homonuclear  diatomic  Lennard-Jones  (LJ)  fluids  for  which  a 

limited  amount  of  Information  from  molecular  dynamics  simulations  Is 

available.  For  two  components  of  approximately  equal  bond  length  but 

different  LJ  parameters  our  results  compare  favorably  with  the, machine 

computations.  From  the  very  few  simulation  data  available  for 

equimolar  mixtures  of  molecules  differing  only  in  their  bond  lengths 

we  cannot  reach  any  firm  conclusion.  Alternative  procedures  for 

treating  general  molecular  fluid  mixtures  are  discussed. 

(Aooeared  tin  Jour.  Chem.  Phys.  ,  15  December  1984)  Reprints  under 
separate  cover  when  available 

8.  New  Systematic  Expansion  of  the  Electric  Field  Distribution 
In  Plasmas 

(Angel  Alastuey,  Carlos  A.  Iglosias,  Joel  L.  Lebowitz  and 
Dominique  Levesque) 


He  derive  a  new  systematic  expansion  of  the  electric  field  distribution 
at  a  test  charge  immersed  in  an  Infinite  two-  or  three-dimensional  one- 
component  plasma.  The  lowest  order  truncation  of  this  expansion  leads 
to  a  mean  field  theory  very  similar  to  the  Adjustable  Parameter 
Exponential  Approximation  (APEX).  The  next  order  corrections  to  this 
mean  field  theory  are  expiicitely  computed  in  terms  of  the  distribution 
functions  of  the  plasma  particles.  All  these  approximations  arc 
compared  to  the  Monte  Carlo  results  for  a  two-dimensional  system  at 
r  "2  and  various  test  charges.  The  systematic  approximations  appear 
to  be  useful.  Even  the  zeroth-order  approximation  is  quite  accurate  fot 
large  test  charges  or  strongly  coupled  systems  and  the  next  order 
Improves  on  it.  still  APEX  is  found  to  be  most  reliable (as  it  is  also 
in  three  dimensions)  and  remains  accurate  in  the  practically  interes¬ 
ting  limit  when  the  test  charge  vanishes*  i.e.  at  a  neutral  atom. 

0 

(Ap^ared  In  Phys.  Rev.  A,  Vol.  30,  No.  6,  November  1984)  Reprints 
under  separate  cover 
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9.  Tho  Two-Dlmonslonal  Ono-Componont  Plasma  In  an  Inhomogonoous 
Dackground:  Exact  Results 

(A.  Alastuoy  and  J.  L.  Lsbowlts) 

Wo  study  the  general  inhomogeneous  two-dlmonslonal  jo  I  Hum  wliorc 
the  background  density  varies  in  one  space  direction  only.  At  r  -  2, 
explicit  functional  representations  of  the  one-  and  two-body  dons  I  tic 
of  the  particles  arc  derived  in  terms  of  tho  electrostatic  potential 
created  by  the  background.  The  present  model  can  be  used  for 
describing  a  large  variety  of  charged  Interfaces. 

CAppeared  In  J.  Physique  45,  pp.  1859-1874,  December  1984)  Roprlntr. 
under  separate  cover 

10.  Van  der  Vaals  One-Fluid  Theory:  Justification  and  Generalisation 
(David  MacGowan,  Joel  L.  Lebowitz  and  Eduardo  M.  Valsmnn) 

We  describe  an  approach  to  Van  dor  Waals  one- fluid  theory 
based  on  thermodynamic  consistency  and  propose  a  method  for  general¬ 
ising  it  to  non-conformal  fluids. 

(To  appear  in  Chem.  Phys.  Letts.)  Preprint  under  separate  cover 

11.  Long  Wave  Length  Oscillations  in  an  Inhomogeneous  One  Component 
Plasma 

(Joel  L.  Lebowitz  and  Ph.  A.  Martin) 

The  perfect  screening  of  charge  fluctuations  in  an  equilibrium 
plasma  is  extended  to  the  time  displaced  structure  function  of 
a  general  inhomogeneous  one  component  plasma.  Wc  find  that  the 
long  wavelength  modes  oscillate  undamped  with  a  single  frequen¬ 
cy  cZ  ,  being  an  angular  average  of  squares  of  plasma  fre¬ 
quencies  cJtm  in  uniform  systems  with  density 

Our  resultl  are  derived  rigorously  from  the  DBGKY  hierarchy  un¬ 
der  some  reasonable  assumptions  on  the  spatial  decay  of  corre¬ 
lations  and  contains  as  special  cases  previously  obtained  results 
of  this  kind. 

Preprint  under  separate  cover 


12.  Angular  Correlations  in  Dense  Hot  Diatomic  Fluids 

(David  MacGowan I  James  Daniel  Johnson  and  Milton  Samuel 
Shaw) 

Molecular  dynamics  (MOI  simulation  data  for  rigid  diatomic  models  of  and  CO^ 
under  conditions  of  extremely  high  density  and  temperature  are  analyzed  for  static 
correlation  functions.  The  results  show  some  significant  qualitative  differences  from  those 
for  diatomic  fluids  at  normal  densities  and  temperatures  (i.e.,  near  the  triple  point).  For  a 
single  thermodynamic  state  of  N^.  the  radial  distribution  functions  (ROFs)  of  the  (spherical) 
RAM  and  median  potentials  are  found,  also  by  MO.  Whereas  the  median  gives  good 
thermodynamic  results  and  poor  centers  correlation  functions.  RAM  produces  just  the 
opposite.  Thus  no  explanation  in  terms  of  distribution  functions  is  found  for  the  success 
of  the  median  for  thermodynamics  although  an  empirical  correlation  is  found  between  the 
breakdown  of  median  thermodynamics  for  CO^  and  a  distinctive  feature  of  the  molecular 
correlation  functions. 

(To  appear  in  Jour.  Chem.  Phys.)  Reprints  under  separate  cover  when 
available 

13.  Exact  Results  for  the  TVfO'Dinwnsianal  Che -(component  Plasma 
(D.  Nicolaides) 

We  show  that  the  free  energy  and  correlation  functions  of  the  2D  OCP  at  the 
special  temperature  T^,  ■  can  be  explicitly  computed  by  solving  a  free- 

field  Dirac  equation. 

(Appeared  in  Phys.  Letts.,  Vol.  103A,  No.  5,  9  July  1984)  Reprints  under  separate 
cover 

14.  Exact  Results  for  the  1Vra*Diffiensional  IVro-Ckxnponent  Plasma 
(D.  Nicolaides) 

Ihe  free  mergy  and  correlation  functions  of  a  two-dimensional  two-component 

Plasna  are  explicitly  conqxited  at  the  special  temperature  >  e^/21cg. 

(Appeared  in  Phys.  Letts.,  Vol.  103A,  No.  1,2,  18  Jime  1984)  Reprints  tnder 
separate  cover 


IS.  Equivalent  Potontlala  for  Equations  of  State  for  Fluids  of  Non- 
spherlcal  Molecules 

(G.  0.  Vllllaina,  J.  L.  Lebo«rltz  and  J.  K.  Percus) 

We  analyze  the  extent  to  which  the  equation  of  state  and  other 
thermodynamic  properties'  of  systems  of  hard  nonsphorlcal  molecules  can 
be  obtained  from  a  density  Independent  hard  sphere  reference  system. 

We  conclude  that  the  median  and  Darke r-Henderson  prescriptions  effec¬ 
tively  reproduce  all  data  now  available.  We  discuss  the  motivation 
for  these  two  formulations  in  detail. 

(Appeared  in  J.  Chem.  Phys.  ,  ^(4),  15  Auftust  1984)  Reprints  under 
separate  cover 

16.  The  Ideally  Polarizable  Interface:  Integral  Equations 
(M.  L.  Rosinberg,  L.  Blum  and  J.  L.  Lebowitz) 

The  integral  equations  used  in  the  microscopic  theory 
of  the  electric  double  layer  arc  extended  to  the  case  of  an 
Impermeable  interface  separating  two  conducting  media.' 

(Ionic  solutions,  or  plasmas) .This  system  is  a  model  for  an  Ideally 
polarizable  inter  face.  Exact  relations  are  given  for  the 
contact  values  of  the  one  particle  density  function, and  also 
for  the  pair  correlation  functions. 

We  solve  numerically  the  Poisson-Bol tzmann  (PB) , 
the  (iypernetted  Chain  (HNC)  and  Hean  Spherical  (MSA) 
approximations, and  compare  the  results  to  the  exact  solution 
of  the  one  component  plasma  in  two  dimensions. 

Preprint  under  separate  cover 
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